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Abstract 



Jets physics in heavy ion reactions is an important new area of active research at the Relativistic Heavy Ion Collider (RHIC) and 



at the Large Hadron Collider (LHC) that paves the way for novel tests of QCD multi-parton dynamics in dense nuclear matter. At 
present, perturbative QCD calculations of hard probes in elementary nucleon-nucleon reactions can be consistently combined with 
P\J the effects of the nuclear medium up to 0{a]). While such accuracy is desirable but not necessary for leading particle tomography, 



OO 



Jt is absolutely essential for the new jet observables. With this motivation, we present first results and predictions to 0{al) for 
the recent LHC lead-lead (Pb+Pb) run at a center-of-mass energy of 2.76 TeV per nucleon-nucleon pair Specifically, we focus 
'on the suppression of the single and double inclusive jet cross sections. Our analysis includes not only final-state inelastic parton 
interactions in the QGP, but also initial-state cold nuclear matter effects and an estimate of the non-perturbative hadronization 
^corrections. We demonstrate how an enhanced di-jet asymmetry in central Pb+Pb reactions at the LHC, recently measured by the 
(— I ATLAS and CMS experiments, can be derived from these results. We show quantitatively that a fraction of this enhancement may 
O |be related to the ambiguity in the separation between the jet and the soft background medium and/or the diffusion of the parton 
shower energy away from the jet axis through collisional processes. We point to a suite of measurements that can help build a 
'consistent picture of parton shower modification in heavy ion collisions at the LHC. 
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.1. Introduction 

One of the most intensively studied processes in high- 
energy physics is jet production, which can be used to accu- 
rately test the fundamental properties of Quantum Chromo- 
Uynamics (QCD) and explore new physics beyond the Stan- 
dard Model lit]. At the parton level, jet observables, such 
as cross sections, jet shapes or event shapes, can be com- 
puted systematically and reliably order-by-order in perturbative 
Quantum Chromodynamics (pQCD). At the hadron level, non- 
perturbative effects related to hadronization and the underly- 
ing event can be estimated and combined with the pQCD re- 
sults to produce theoretical prediction suitable for comparison 
to the experimental measurements IH 0, Isl Si- The highest 
center-of-mass energy hadronic collisions available to date at 
the CERN Large Hadron Collider (LHC) guarantee an abundant 
yield of large transverse momentum hadrons and jets and pro- 
vide, through experimental measurements and constraints on 
theoretical models l|5n^, an ideal opportunity to advance per- 
turbative QCD. 

In relativistic heavy ion collisions at the LHC, jet and par- 
ticle production will be altered by the collisional and radiative 
final-state interactions of the parent partons which propagate 
through a region of hot and dense deconfined matter, or the 



quark-gluon plasma (QGP). It has been demonstrated that a sig- 
nificant part of the partons' energy can be lost due to medium- 
induced gluon bremsstrahlung [Tj], and uncertainties in energy 
loss calculations have been commented upon in liSfl. The pre- 
dicted related suppression of leading particles at the Relativistic 
Heavy Ion Conider(RHIC) and at the LHC 13) is now well es- 
tablished ifioi [Tl , 121. including its non-trivial momentum de- 
pendence verified by the ALICE Collaboration. On the other 
hand, the modification of jets in dense QCD matter requires 
much deeper understanding of the in-medium parton shower 
production mechanisms lll3ll . For inclusive jets, tree-level cal- 
culations can be used to discuss qualitatively the effects of final- 
state interactions but 0(al) results il4il are necessary for quan- 
titative comparison to the preUminary jet experimental data in 
gold-gold (AuH-Au) and copper-copper (Cu+Cu) reactions at 
RHIC [11 EH- 

For tagged jets and di-jets lowest order, 
for example (9(a^), 0(Gfas) and OiOemO:.,), calculations fail to 
describe even qualitatively the double differential momentum 
distributions llST and one-loop pQCD results are required. 
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Recently the ATLAS and CMS Collaborations reported a 
significant enhancement in the transverse momentum imbal- 
ance of di-jets produced in central Pb+Pb collisions at the 
LHC dlMl- It has been qualitatively argued that such asym- 
metry A J may be a consequence of jet quenching in the QGP 
created in heavy ion reactions at LHC 112 ill . Attempts to explain 
its magnitude based on Monte Carlo simulations with a Pythia 
generated p+p baseline have also been presented l22l l23l 12411 . 
In Ref. 12511 the important observation was made that the huge 
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heavy ion background and its fluctuations may significantly en- 
hance the di-jet asymmetry distribution. 

In this Letter we report first results and predictions at 0(al) 
for inclusive jet and di-jet productions in central Pb+Pb colli- 
sions with y/Sffff = 2.76 TeV at the LHC. We combine consis- 
tently the 0{a]) calculations in nucleon-nucleon reactions, val- 
idated through comparison to the Vs = 7 TeV p+p data at the 
LHC, with parton energy loss eff^ect in the hot QCD medium 
at Oia^^af^) lH]. 0{al), Oia^.af^) results represent next-to- 
leading order accuracy for inclusive jet production but for the 
more differential observables, such as the di-jet asymmetry at 
Ay > 0, this may be the first order that gives a non-trivial (9^0) 
contribution. We demonstrate the relation between the cross 
section attenuation for inclusive jets and di-jets and the result- 
ing enhancement in the di-jet asymmetry A j. We find that a 
part of this enhancement may arise from the uncertainty in the 
jet/background separation and/or the dissipation of the energy 
from the shower into the QGP due to collisional processes 127 1 . 
To further differentiate these two possibilities is subject of fu- 
ture work. 

Our Letter is organized as follows: we present benchmark re- 
sults for inclusive jet and di-jet production in ^fs - 7 TeV p+p 
collisions at the LHC in section |2l Results and predictions for 
the central Pb+Pb reactions at the LHC at yfJim = 2.76 TeV are 
given in section[3] Our conclusions are presented in section |4] 



2. Inclusive jet and di-jet production in p+p collisions 

Within the framework of the collinear perturbative QCD fac- 
torization approach the lowest order (LO) single and double in- 
clusive jet production cross sections can be obtained from: 
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In Eq. ([TJ 5 = {Pg + Phf- is the squared center-of-mass energy 
of the hadronic collision and d^Ej - ErdErdtp with (f being 
the azimuthal angle of the jet. 



are the lightcone momentum fractions of the incoming partons. 
We denote by (pa,b/N{xa,b) the distribution function of partons 
a,b in the nucleon. |M„fc^,yp are the squared matrix elements 
for a + b i + j processes, for massless partons see 112 811 . In 
the double collinear approximation the jets are produced back- 
to-back and identified with the parent partons. Formally, the 
partonic cross section in Eq. ([TJ is convolved with the jet phase 
space S2{p'{,p^',) (dimensionless in our convention). At tree 
level, this phase space factor is simply a sum of products of 
^-functions. For inclusive jet production, the second jet is inte- 
grated out. 
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• ATLAS Inclusive Jet Data 
— NLO pQCD X Non-Pert. Correction 
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Figure 1 : Comparison between the i? = 0.4 inclusive jet cross section measured 
by the ATLAS Collaboration in p+p collisions at = 7 TeV and the 0(a^^) 
QCD theory which includes non-peiturbative conections. Expeiimental en'or 
bars are statistical only. 



The first obvious deficiency of the LO formahsm is its in- 
ability to describe the jet radius dependence of the physical ob- 
servables Il4i1 . The second much more significant shortcom- 
ing is the inadequate description of the differential distribu- 
tion of tagged jets |fl8| and di-jets due to the 5-function con- 
straint in Eq. ([TJ. At 0{al) and beyond jets are defined by 
the jet-finding algorithms characterized by a radius parameter 
R = yjlS.y'^ + A)p~, which gives the jet's spacial extent in rapid- 
ity y and azimuthal angle 0. Generally, the cross sections of the 
jet production at the 0{a]) can be expressed as ll29ll : 



do- _ I r 
dV„ ~ T.J 

4/ 



do-{2 -> 2) 
dy\dy2d^ET\d^ET2 



dVi 



do-{2 3) 



dyidy^dy^d^Ej id^ET2d^ET3 



(2) 



where V,, = dy\ ■ ■ ■ dy„d^ET i ■ ■ ■ d^Ej,, represents the phase 
space for multi-jets or multi-partons, respectively. The first 
term on the right-hand side of Eq. (|2]i gives the contribution 
from 2^2 processes, including the 0{a^^) virtual correc- 
tions. The second term in Eq. (|2]i denotes the contribution from 
2 — » 3 processes. The jet size R and algorithm dependence 
are contained in the function S T,{p^^, p'^, p'i) lITiil . which is sig- 
nificantly more complicated than S 2(1^1, P2)- In the calcula- 
tions that follow we will take advantage of the Ellis, Kunszt 
and Soper (EKS) numerical implementation of 0(al) jets pro- 
duction in hadron-hadron collisions 1I29I1 and the CTEQ6.1M 
parametrization set of parton distribution functions 1139 1. 

Measurements of inclusive jet production in p+p collisions 
at yfs - 7 TeV with jet cone radii R - 0.4, 0.6 and in different 
rapidity regions have been carried out by the ATLAS Collab- 
oration at the LHC Perturbative QCD calculations can be 
corrected to the hadron level using parametrization of universal 
non-perturbative effects fyp{ET,R) extracted from Monte Carlo 
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simulations as follows: 

j^hadi-on 



^^parton 



dEj 1 • • ■ dEj „ dEj 1 ■ ■ ■ dEj „ 



Y[hAETi,Rd- (3) 



In this Letter we employ the parametrization of non- 
peiturbative effects provided by ATLAS jsl in all calculations 
except the ones presented in the bottom panel of Fig. |6] The 
only observables that are noticeably affected by ff^p are the ones 
that involve too significantly different radii Ri, R2- Returning 
to jet production in p+p collisions, we find that an excellent de- 
scription of the differential jet spectrum is obtained in the mea- 
sured region of Ej between 50 GeV and 500 GeV where the 
yield falls by more than 5 orders of magnitude. We present one 
example for rapidity \y\ < 0.3 and 7? = 0.4 in Fig. [1] Variation of 
the factorization and renormalization scale - - Ej in the 
interval (Et/2,2Et) has less than 10% effect on the inclusive 
jet cross section. 

The evaluation of di-jet production in hadronic reactions is 
slightly more involved. At 0{al) the naive exact transverse en- 
ergy equality Eji - Eji between the two jets is broken by 
2^3 parton splitting processes. Loop corrections in 2 — > 2 
processes, proper jets definitions and finite Ej bins ensure well- 
behaved physical cross sections. To quantify the resulting im- 
balance we modify the 0{al) EKS code 12911 and define the 
dimensionless: 
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for brevity. In Fig. |2]we show the three-dimensional (3D) plot 
of (T as a function of the transverse energy of the two jets Ej i 
and Et2 in the rapidity region \y\ < 2.8 for the same jet size 
R - 0.4. Transverse energy bins AEj \ , AEt2 - 10 GeV have 
been chosen. While the cross section reaches its maximum for 
Eti ~ Et2, its most striking feature is how broad and slowly 
varying the di-jet yield is away from the main diagonal. As 
we will see below, such broad distribution limits the additional 
asymmetry that jet propagation in the QGP will induce. 

One particular distribution that can be extracted from our 
general di-jet cross section calculation is that of the di-jet asym- 
metry A J, defined as 



A, = 



Et 1 — E 



T2 



Ej \ + E 



(5) 
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Making a change of variables from {Ej i , £7 2) to (Ay, Ej 2) with 
Eji - £'7-2(l+Ay)/(l-Ay) and then integrating over3'i,>'2,£r2 
we can express the Aj distribution from the di-jet Eji, Ej2 
spectrum as follows: 



dAj 



dE 



T2 



2Et2 do-[ETi{Aj,ET2)] 



(6) 



(I - Aj)^ dy\dy2dET \dET2 

where £7 2 max is given by the jet production kinematics. By 
evaluating Aj > we ensure that £71 > Et2- Note that at LO 
we always have Ay = 0. 
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Figure 2: A three-dimensional representation of the Oia^) di-jet cross section 
(T, defined in Eq. ^4), in = 7 TeV p+p collisions at the LHC versus Ej 1 and 
£7-2. We have used a jet radius R = 0.4 and a rapidity interval |v| < 2.8. 
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Figure 3: The noiTnalized differential distribution of the di-jet asymmetry A j, 
defined in Eq. (3), is evaluated at Oia],) and compared to the ATLAS experi- 
mental measurement at = 7 TeV. Solid (dashed) lines include (do not in- 
clude) non-perturbative corrections, respectively. 



With the 0{al) numerical results for & at hand, we can easily 
compute cr^^dcr/dA j from Eq. (|6]l and confront it with the ex- 
perimental data. The ATLAS Collaboration has measured the 
di-jet asymmetry distribution in p+p reactions at = 7 TeV 
with the following kinematics cuts: Eji > 100 GeV for the 
leading jet and £'y2 > 25 GeV for the subleading jet ||T§tl. The 
corresponding comparison to the perturbative QCD theory is 
presented in Fig. [3] and shows excellent agreement. Variation 
of the factorization and renormalization scale fj.f - fiR - Ej in 
the interval {Et/2, 2Et) has less than 2% effect on the A y dis- 
tribution. The dashed line represents the parton level 0{a].) re- 
sults and the solid line includes non-perturbative effects, which 
are seen to be practically negligible for this observable. Strict 
transverse energy ordering Et2 < Et\ for the leading and sub- 
leading jets implies dcr/dA y = when Ay = 0. 
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Fixed 0{a^^) calculations of dcr/dAj will diverge for Ay — » 0. 
To ensure proper cancellation of real and virtual contributions 
to the physical cross section, these evaluations rely on finite bin 
sizes in transverse energy, rapidity and azimuth. For example, 
in numerical simulations based on the EKS code, results are 
calculated at random positions in phase space and any required 
grid for the differential cross sections is given by interpolation 
with a Gaussian smearing function ll29ll . Thus, our numerical 
simulations yield continuously smooth results of Eji ~ Eti 
or, equivalently, Aj ^ Q. On the other hand, fine resolution 
cannot be achieved when £7 1 Etioi Aj ^ Q. To investigate 
this region, resummation beyond fixed order calculations will 
be needed. Q 



3. Inclusive jet and di-jet production in heavy ion collisions 

The bases for the evaluation of multi-jet cross sections in 
heavy ion collisions are the corresponding cross sections in the 
more elementary nucleon-nucleon reactions. In Ref. llsoll it was 
shown that in QCD the final-state process-dependent medium- 
induced radiative corrections factorize in the production cross 
sections for the example of a single jet. The generalization to 
multiple jets proceeds as follows 



do-{ei , 



«-jet 
quench. 



dcr{€i, 
®P„{e„) \Ji{ei)\ 



5Pi(ei)- 
|/„(e„)l . 



(7) 



Here, P,(e,) is the probability that the /th jet will lose a frac- 
tion e, of its energy, |y,(e,)| are phase space Jacobians, and ® 
denotes an integral convolution. Any dependence on the jet re- 
construction parameters is not shown explicitly in Eq. (|7). Fur- 
thermore, possible small contributions of the medium-induced 
parton shower associated with one jet to other jets that are very 
close in phase space (y, (p) are not included. For inclusive and 
approximately back-to-back jets that we consider in this Letter, 
this effect is negligible. 

The application of Eq. (|7]i to inclusive jet production has been 
discussed extensively in Refs. iflslfTill . where the angular dis- 
tribution of the medium-induced gluon radiation is computed 
using the Reaction operator approach J?! IssJ. Here, we sum- 
marize the main result. The medium-modified jet cross section 
per binary scattering is evaluated as 



1 dcr^{R) 



{Nbm) dyidEji " Jfi=o 



dei PnJeuETi) 



, CNM.NLO, c., ^ 



In Eq. ([8]l the phase space Jacobian reads 

l^,(e,)| = l/(l-[l-/W-,prMf')> 



(8) 



and E'j. - \ Ji{ei)\ETi in the argument of cr™'^-'^^°. Note that in 
the expression above we have denoted by 



du) 



dcodr 



lo' drj^'^'^doj 



dojdr 



(9) 



the part of the fractional energy loss e, that is redistributed in- 
side the jet. Thus, /(/?,, Pjf)q,g plays a critical role in defin- 
ing the contribution (or lack thereof if /(/?,, Pjf)q,g — > 0) of 
the medium-induced bremsstrahlung to the jet. Note that in 
Eq. (|9j Pjf is a parameter that can be used to simulate pro- 
cesses that can alter the energy of the jet beyond medium- 
induced bremsstrahlung. For example, by choosing pjf > 
GeV one can investigate phenomenologically the uncertainty 
in the jet/background separation and the diffusion of the parton 
shower energy away from the jet axis due to collisional pro- 
cesses. Specifically, one can adjust the value of p™" to reflect 
the collisional energy loss of the parton shower {6Et) evaluated 
in Ref. H. 

Cold nuclear matter (CNM) effects prior to the QGP forma- 
tion, such as initial-state energy loss, coherent power correc- 
tions and the Cronin eff'ect Bill , can be incorporated in the cross 

sections for quark and gluon jet production '^'^.^^^ and their 
multi-jet generalization. At the LHC we are interested in jets 
of Ej > 20 GeV and at these energy scales of the efl'ects listed 
above only initial-state energy loss, which is also derived in 
the framework of the Reaction operator approach, may play a 
role II32I I33II . In this work we evaluate the effects of cold nu- 
clear matter on jet production in Pb+Pb collisions at the LHC 
and comment on their relevance to specific observables. 

For di-jets produced in heavy ion reactions, the correspond- 
ing cross section can be expressed as 



1 



do^'^iR) 



{Nbin > dy idy2dET\dET2 



V I c/ei I de2 



qg.ig-gg 



(1 - [1 - f{Rupf^)q,g\ei) (1 - [1 - /(-R2,p?'^),,,]e2) 

(10) 



, CNM,NLO,^, ^, . 
qq,qg,gg \^Tl' T2' 



' We thank the reviewer for a suggestion to clarify this point. 



dyidy2dE^ ^dEj2 

The evaluation of the inclusive and di-jet cross sections at 
the LHC proceeds as follows: the hard jet production processes 
are distributed in the transverse plane according to the binary 
collision density o-i„(dTAAib)/dx±). The density of the QGP is 
taken to be proportional to the participant density dNpait./dx±. 
Through isospin symmetry, necessary to account for the neutral 
particles, and parton-hadron duality lllSll the gluon rapidity den- 
sity in central Pb-i-Pb collisions at the LHC at -sfsNN - 2.76 TeV 
is estimated to be dN^/dy - 2650 from the ALICE experimen- 
tal measurement of dN'^'^drj fsill . We take into account the lon- 
gitudinal Bjorken expansion of the medium and, by assuming 
local thermal equilibrium, we can evaluate the physical quanti- 
ties, such as the Debye screening scale, the in-medium jet scat- 
tering cross sections, and the mean free paths [Tsl 14 1, that are 
employed in the calculation of the radiative energy loss ll32ll . 
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The energy loss calculation, which uses the collision geom- 
etry described above, includes kinematic constraints, and pro- 
vides the fully differential distribution of the medium-induced 
bremsstrahlung ll32ll . is computationally demanding and per- 
formed separately. Using the analytic leading power depen- 
dence (up to logarithmic corrections) of the energy loss on the 
coupling between the jet and the medium ^n,ed and the path of 
the parton propagation in the QGP: 



4 

?med 



f 

J to 



fp(xo± + nj_f, f) dt , 



(11) 



we can precisely account for/study the sensitivity to the de- 
tails of the in-medium jet dynamics by scaling the averaged 
bremsstrahlung spectra. Note that g^^^/'^^ above should not 
be confused with a'^ in the medium-induced bremsstrahlung 
vertex. In Eq. ( fTTT i xq^ is the jet or di-jet production point 
in the transverse plane and is the direction of its propaga- 
tion. Finally, we evaluate in real time the probability distri- 
bution for the parent partons to lose fractions of their energy 
= Zin(w„)/7£/ due to multiple gluon emissions in the Poisson 
approximation tgiJ; 



PaJeddei = 1 



Jo 



eiPqAeddei 



AEi 



(12) 



This approximation is consistent with the inclusive gluon emis- 
sion calculations currently used in jet quenching phenomenol- 
ogy IBS. 

3.1. Quenching of inclusive jets 

To investigate jet production in relativistic heavy ion colli- 
sions and quantify its deviation from the baseline results in el- 
ementary hadron-hadron reactions we define a generalized nu- 
clear modification factor for multi-jet final states as follows: 

rin-jet^P jr . p p „inin „™n\ 



dcr^ (ETi-ET„;R\-R,„PT'r- -p't'" ) 



dcr"i'(ETvET„;Ri-R:„p^';'-pT:;) 
dyydy'i^dEj \ ■■■dEj^, 



(13) 



Next, we calculate numerically the cross sections and nuclear 
modification factors E^^^^^ for inclusive jet and di-jet production 
(« = 1,2) in central Pb+Pb collisions at LHC with yfsNN - 
2.76 TeV. 

Fig. |4] illustrates the dependence of the suppression ratio 
Rjf^^^^ for inclusive jets in central Pb+Pb collisions on the jet 
size R. The top panel shows the net suppression resulting from 
both initial-state and final-state effects for R - 0.2, 0.6 and the 
bands represent the variation in coupling strength between the 
jet and the medium. The dashed, solid, and dot-dashed lines 
correspond to ^^ed =1.8, 2, 2.2, respectively. Note that for 
small radii the suppression of jets can be comparable in magni- 
tude to the suppression of inclusive particles. This is illustrated 
by including a reference charged hadron quenching measure- 
ment by the ALICE Collaboration lll2ll in Fig. |4] Preliminary 
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Figure 4: Top panel: the Ej dependence of the nuclear modification factor for 
different jet cone sizes i? = 0, 2, 0.6 is calculated in central Pb+Pb collisions 
at the LHC -\jSNN = 2.76 TeV. Bands represent the variation in the couphng 
strength between the jet and the medium. Bottom panel: the relative contribu- 
tion of cold nuclear matter effects to is illustrated for i? = 0.2, 0.4, 0.6, 0.8. 
ALICE experimental data on chai'ged hadron suppression in central Pb+Pb col- 
lision is shown for reference. Prehminary ATLAS Rcp data is also shown for 
both R = 0.2, 0.4. 



ATLAS data 11361] on i^^p*^', a ratio similar to R]^^^'^ where pe- 
ripheral A+A reactions are used as a baseline, has now been in- 
cluded in the figure^. The bottom panel of Fig.|4]presents a cal- 
culation for R = 0.2, 0.4, 0.6, 0.8 and for a fixed gmed =2 with 
(solid line) or without (dashed line) cold nuclear matter effects. 
For fixed centrality, CNM effects, here dominated by initial- 
state energy loss, do not depend on the jet size or jet finding 
algorithm and become more relevant, relatively speaking, for 
large radii R. Even though on an absolute scale this additional 
suppression is not very large, it is more significant in compar- 
ison to the or Dell-Yan production processes 1 18 371. 
These latter processes are dominated by ^ + ^ initial states 
and jet production discussed in this manuscript arises primar- 
ily from g + g (and g + q{q) at larger Ej) processes. 

Initial-state CNM effects in heavy ion collisions can be min- 
imized by taking the ratio of jet cross section at two different 
radii [do-{Ri)/dET]/[do-(R2)/dET] Q. Since the size R de- 
termines what fraction of the parton shower is reconstructed 
as a jet, it affects the jet cross section. In heavy ion reactions 



-Two weeks after this manuscript was been posted to the archive, the AT- 
LAS Collaboration showed preliminary results on inclusive jet quenching in 
central Pb+Pb collisions at the Quarii Matter 2011 conference. We include 
two data sets to give the reader a feel as to how the predictions presented here 
compare to the new data. 
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Figure 5: Ratio of the inclusive jet cross sections in central Pb+Pb collisions 
at LHC at = 2.76 TeV for two dilferent radii R, =0.2 and R2 = 0.4. The 

bands show results with different extrapolation of non-perturbative corrections 
to small radii. The lines show effect of different coupling strength between the 
jet and the medium. 



the cone size dependence is amplified by the fact that medium- 
induced parton showers have a broad angular distribution in 
comparison to the ones in the vacuum yjl]- This is shown in 
Fig. |5]for Ri - 0.2, R2 - 0.4 and the dashed, solid, and dot- 
dashed lines correspond to three different ^,ned = 1-8, 2, 2.2. 
As the radius varies, specific non-perturbative effects, unfortu- 
nately, become more important. Typically, they are expressed 
as an average momentum shift |0, [sstl ™d related to "splash- 
out" hadronization effects and "splash-in" initial-state radia- 
tion/background contribution: {dpr) = A/R + BR^. The phys- 
ical effect of a momentum shift is to alter the measured cross 
section and this change can be isolated in a multiplicative fac- 
tor list]. Since background effects are the dominant uncer- 
tainty in jet physics with heavy ions, we will discuss them 
separately. With this motivation, we consider a hadronization- 
motivated extrapolation of the ATLAS parametrization of non- 
perturbative effects to small radii: fyp - a + b/R. The ap- 
plication of this non-perturbative correction to the calculation 
of [do-(Ri = 0.2) /dET]/[do-(R2 = 0.4)/dET] in central Pb-hPb 
collisions at the LHC is shown by the cyan band in Fig.|5] where 
we assume that for jet production at large Ej the "splash-out" 
hadronization effects in A+A collisions are the same as that 
in elementary p-i-p collisions with the understanding that the 
hadronization time for jets at high Ej is longer than the dura- 
tion time of the QGP. Note that the non-perturbative effect can 
change significantly the cross section ratio relative to the 0(a]) 
parton level result for small R. It is, therefore, critical to con- 
strain its magnitude as accurately as possible in the simpler p+p 
reactions. Moreover, higher order corrections may affect these 
ratios to some extent and a recent study has shown that in p-i-p 
collisions higher order correction may lower the jet cross sec- 
tion ratios considerably 113 811 . 

Preliminary RHIC results suggest that the jet size depen- 
dence of jet attenuation may have already been observed in 
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before we discuss di-jet production in heavy ion reactions, we 



Figure 6: The effects of parton shower collisional energy loss or misinter- 
preting some of the jet's energy as soft background on R^'^' (top panel) and 
<t(R[)/<t(R2) (bottom panel) are simulated with p™" = 20 GeV and repre- 
sented by a cyan band. Results are contrasted to the p™'" = GeV case, 
represented by a yellow band, which is no collisional energy loss or unam- 
biguous jet/background separation. Preliminary ATLAS Rcp data is shown for 
R = 0.2, 0.4. 



comment on the difficulties related to the measurement of jet 
observables. In central Pb-i-Pb collisions at the LHC, for a typ- 
ical jet size R = 0.5, on the order of 100 GeV of the jet energy 
is interpreted as QGP background and subtracted from the total 
reconstructed energy ll20ll . While a simple jet+uniform back- 
ground model appears reasonable in heavy ion reactions, it is 
not based on first-principles theory. In what follows we demon- 
strate the consequences of misinterpreting less than 20 GeV of 
the jet energy redistributed by the QGP medium inside the jet 
as uncorrected soft background. This is < 20% of the typi- 
cal subtracted Er and in our approach lll3ll can be simulated 
by choosing p™'" = 20 GeV in Eq. (|9|i- We denote the actual 
amount of energy by {SEj) and emphasize that it depends on 
the jet radius R. More importantly, we recall that a recent cal- 
culation of the energy transmitted by a parton shower to the 
medium |l27tl A^Cshower ^ QGP) found that for LHC con- 
ditions A£'(shower — > QGP) - 20 GeV is well within reach, 
especially for a gluon-initiated shower While the growth in 
the rate of collisional energy loss is due to the proliferation of 
medium-induced partons, the rate for each individual parton is 
calculated to Oig^^^^. Of course, one should keep in mind that 
experimental observables are sensitive to the physics effects on 
the fraction of the in-medium parton shower that falls within the 
cone radius R. We first choose p™" - 20 GeV to simulate the 
collisional energy loss effects evaluated in |27] on the parton 
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shower for a very large radius R ^ 2. We then vary R to obtain 
the coUisional energy loss effect on the part of the parton shower 
that constitutes the jet. For example, for a 100 GeV gluon jet 
choosing p™" = 20 GeV in Eq. ® leads to {6Et} ^ 15 GeV 
for R = 0.5 and to {6Et} ^ 5 GeV for R = 0.2. The effect 
of collisional energy loss decreases quickly with the jet radius 
since the medium-induced parton shower has a characteristic 
large-angle 13511 distribution and most of it is already outside of 
the jet cone for small R |[l3l [l4l[l8l. 

The result of our simulations is shown in Fig. |6]for gmed = 2, 
where the default choice p™" - GeV is illustrated by a yellow 
band and the choice p™" = 20 GeV - by a cyan band. In the top 
panel, the clear dependence of R^/^'^ on the jet size, exemplified 
hy R = 0.2, 0.4, 0.6, practically disappears if a sizable fraction 
of the particles created by the medium-induced bremsstrahlung 
diffuse outside of the jet due to collisional processes or are sub- 
tracted as soft background. It can easily be seen that in this 
case and in the absence of strong /^-dependent non-perturbative 
corrections cr{R\) j a-{R2) ~ 1. This is demonstrated in the bot- 
tom panel of Fig. |6]for three gmed = 1-8, 2, 2.2 and the result 
is very different in comparison to letting = GeV. We do 
not include the non-perturbative correction in the bottom panel 
of Fig. |6] since it can result in noticeably different cross sec- 
tions for significantly different /?i , Ro and will distract from the 
difference between the simulations that only include radiative 
energy loss processes and the ones that include the effects of 
the collisional interactions of the parton shower We conclude 
that the ambiguity in the jet/background separation and the pos- 
sibility that some of the parton shower energy can diffuse away 
from the jet axis at angles greater than R through collisional 
processes 112711 in heavy ion reactions can significantly affect 
the experimental jet cross sections. 

3.2. Quenching of di-jets and their Ej asymmetry 

We now present our results for the modification of di-jet pro- 
duction in central Pb+Pb collisions at ^jSf4N - 2.76 TeV The 
cross section is evaluated using Eq. (fTOl i. First, we give a 3D 
representation for the suppression ratio R'^^'^^ of the di-jet cross 
section at 0{al) in Pb+Pb collisions relative to the theoretically 
computed pH-p baseline in Fig. |2l In this example we use jets 
with identical cone size R - 0.2, a coupling strength between 
the jets and the medium ^med = 2, and no transverse energy 
cut. The suppression involving both cold nuclear matter ef- 
fects and final-state quenching effects in the QGP is the largest 
along the diagonal Et\ = Et2- In the regime Eji » Eji 
01 En <*; Et2 there is a striking enhancement that can reach 
values much larger than 4 (the maximum value shown no the z 
axis in Fig. [T]). Such enhancement was first identified theoret- 
ically on the example of tagged jets ifisll and is easy to un- 
derstand. The final-state inelastic interactions of jets produced 
in heavy ion collisions redistribute them to smaller transverse 
energies and these jets fill in a region of phase space where the 
cross section is steeply falling, see Fig.|2] The second important 
observation is how broad the region of approximately constant 
TJ^J*^' in the {Et\,Et2) plane is before it hits the strong di-jet 
enhancement. This behavior limits the features in the heavy 
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Figure 7: 0(a^^) calculation of the di-jet cross section suppression ratio 
versus Ej i and Ej 2 in y-s-jvyv = 2.76 TeV central Pb+Pb collisions at the LHC. 
Our result is for Ri = R2 = 0.2 and includes cold nuclear matter effects. 



ion A J enhancement that can be attributed to the physics of jet 
quenching. 

Next, the di-jet asymmetry is computed using Eq. (|6]l and 
compared to recent experimental results published by AT- 
LAS {H and CMS ^ in Fig. [8] and Fig. |9l The ATLAS 
and CMS experiments apply different transverse energy cuts 
for their leading and subleading jets: (Ej 1 > 100 GeV, Et 1 > 
Et2>25 GeV) and {En > 120 GeV, Eti > Et 2 > 50 GeV), 
respectively. Two different calculations are performed and 
shown in the top and bottom panels of Figs.[8]and|9] In all cases 
the baseline p+p asymmetry at = 2.76 TeV is shown by a 
sohd black line. Note that we present the normalized asymme- 
try result (1 /cr)dcr/dA j. 

Fig. [8] shows the sensitivity of the di-jet asymmetry to the 
coupling strength g^ed of the jets to the QGP medium. We use 
one jet size R = 0.4 (we will vary R later). This = GeV 
case is represented by green curves. Qualitatively, we find a 
broadening of the di-jet asymmetry distribution but the depen- 
dence of this broadening to the variation of ^med = 1-8, 2, 2.2 
is quite modest. Quantitatively, approximately 1/3 to 1/2 of the 
extra broadening that experiments attribute to jet interactions in 
matter can be described. Next, we consider colhsional interac- 
tions of the parton shower or an ambiguity in the separation of 
the jets from the background in central Pb+Pb reactions at the 
LHC. This is represented and modeled by p™" = 20 GeV and 
shown by red lines in Fig. |8] In this case, a significantly larger 
broadening is observed. However, only part of it is due to radia- 
tive jet quenching processes. The remainder may be due to re- 
moving the soft particles originating from the medium-induced 
parton shower inside the jet as uncorrected soft background. 
In this sense, our result is similar to the argument presented 
in II25II that background fluctuations can generate much of the 
asymmetry observed by ATLAS and CMS. Alternatively, there 
is the possibility that this energy is dissipated outside the cone 
through collisional processes f27\, something that we will in- 
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Figure 8: Di-jet asymmetry distributions for different coupling strength! in cen- 
tral Pb+Pb collisions at = 2.76 TeV compared to data. The top panel is for 
leading jets of En > lOOGeV, subleadingjets of £7-2 > 25 GeV (ATLAS (H], 
R = 0.4). The bottom panel is for leading jets of £"7- 1 > 120 GeV, subleading 
jet of £72 > 50 GeV (CMS gl. R = 0.5). Black lines are the results for 
p+p collisions under pQCD theory. Green hnes assume perfect jet/background 
separation and are denoted p™'" = GeV. Red lines, denoted p'™" = 20 GeV, 
show the consequences of the ambiguity in the jet/background separation or the 
diffusion of the in-medium parton shower energy away from the jet axis due to 
collisional processes. 



vestigate further in the future. This scenario is also modeled 
by - 20 GeV and in this case we have made quantitative 
connection to theoretical simulations ll27ll . We finally point out 
that even if all the energy associated with the medium-induced 
parton shower is removed, the resulting A j distribution is flat. 
Specifically, a peak in this distribution at finite Ay = 0.3 - 0.4 is 
not compatible with realistic jet quenching calculation^ This 
is easy to be understood from the 0{al) results presented here. 
The interested reader can analyze Fig. |7]and see that the flat A j 
distribution is related to the very broad approximately constant 

p2~jet 

Fig. |9] shows the dependence of the di-jet asymmetry on the 
jet radius R and the the momentum cut p 
pling strength ^^ed = 2. Note that for p 
a significant dependence on the jet size. For pj 
the dependence on the radius is practically eliminated. This 
observation is compatible with the comment by the ATLAS 
Collaboration lfl9ll that their asymmetry measurement has lit- 
tle sensitivity to the choice of R over a wide range of cone sizes 
R = 0.2 - 0.6. It once again stresses the importance of consider- 
ing the interaction between the full parton shower and the QGP 



™" with a fixed cou- 
™" = GeV there is 
™" - 20 GeV, 



^At the recent QM2011 conference, after the results presented in this 
manuscript were posted to the archive, the ATLAS Collaboration presented a 
re-analysis of their A j measurement with improved background subtraction. 
The peak at Ay =: 0.4 is no longer present and the new distribution is consistent 
with flat up to the kinematic bound. 
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Figure 9: Same as Fig [8] but for a fixed jet-medium interaction strength 
gmcd = 2. Dot-dashed, solid, and dashed hnes correspond to S = 0.6, 0.4, 0.2 
respectively. Note that if some of the subtracted "background" is in fact re- 
lated to the medium-induced parton shower or energy diffuses away from the 
jet axis due to collisional processes, the enhanced asymmetry becomes radius- 
independent. 



beyond the medium-induced parton splitting processes 1127 1 
and/or the ambiguity in the experimental jet/background sep- 
aration 10]. 

4. Summary and Conclusions 



Jet production and modification in high-energy nucleus- 
nucleus collisions have been proposed as new ways to unravel 
the properties of the hot and dense QCD medium and to elu- 
cidate the mechanisms of in-medium parton shower forma- 
tion Il3lll4ill8ll . Measurements of jets in heavy ion collisions 
have now become available \vh . With this in mind, in 

this Letter we presented first 0{al), (9(a^a"'') results for the 
single and double inclusive jets production rates in p-i-p and 
central Pb-+Pb collisions at the LHC at center-of-mass energies 
per nucleon-nucleon pair yfsNN - 7 TeV and 2.76 TeV, respec- 
tively. We placed emphasis on the good agreement between the 
perturbative QCD theory and the experimental measurements 
in p-Hp collisions. The importance of reliable 0{al) calculation 
cannot be under emphasized, since for tagged jets and di-jets an 
inaccurate description of the p+p baseline at LO can produce 
even qualitatively incorrect jet quenching predictions lIlSll . 

We first addressed the suppression of inclusive jet produc- 
tion in central Pb+Pb collisions at the LHC and found that it 
is dominated by final-state inelastic parton interactions in the 

32 . 33ll. even though 



QGP. Cold nuclear matter effects 

larger than those for inclusive TP production II37I1 . still do not 



contribute significantly to the attenuation of jets for R < 0.4. 
On the other hand, there is a clear dependence of R^l'^^ on the 
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coupling strength between the jet and the medium gmed- We 
found that for small radii the suppression of jets matches onto 
the suppression of leading particles lll4ll . here represented by 
the ALICE charged hadron measurement. Cold nuclear matter 
effects can be practically eliminated by taking the ratio of the 
cross sections at two different jet sizes o-{Ri)/cr{R2), but in this 
case there is a sensitivity to the non-perturbative hadronization 
effects. We used an extrapolation of the ATLAS parametriza- 
tion of non-perturbative effects to show that for small radii 
R ~ 0.2 they can alter the cross section ratios by as much as 
50%. It is important in the future to extend the experimental 
measurements to lower Ej and better constrain the magnitude 
of the non-perturbative effects. 

Theoretically a clean calculation that focuses on the medium 
modification of jets is easy to perform. Experimentally, the 
situation is much more complicated in heavy ion reactions 
where there is an enormous soft background that may contribute 
~ 100 GeV to the Et of an 7? = 0.5 jet at the LHC. The problem 
of this background separation is both technical and conceptual 
since there is no unified first-principle understanding of heavy 
ion dynamics at all momentum scales. In our calculation we 
investigated the experimental consequences of this ambiguity 
by simulating a subtraction of a radius-dependent {5Ej) from 
the medium-induced parton shower as a part of the soft uncor- 
rected background. Such amount of energy {SEj) may also 
easily be transferred from a 100 GeV shower to the medium 
through collisional processes. The first possibility emulates en- 
ergy loss. The second possibility is the collisional energy loss 
of the parton shower that falls inside the jet cone radius R. We 
used the calculation in Ref 12711 to constrain its magnitude. We 
found that such subtraction wipes out the jet size dependence 
of the quenching observables. 

We also presented the first calculation of the di-jet modifi- 
cation ^^J*^' in heavy ion reaction and found that it is char- 
acterized by a broad suppression region near £71 - Eji and 
strong enhancement for Ejx » £7-2 or £71 <^ Eji- The re- 
sulting 0{al) di-jet cross section was used to evaluate the di-jet 
asymmetry A j. Assuming perfect jet/background separation we 
found that less than 1/2 of the broadening can be accounted for 
the radiative jet quenching calculation. We found little sensitiv- 
ity of the A J observable to the strength of the coupling between 
the jet and the medium and practically no dependence on the 
non-perturbative hadronization corrections. On the other hand, 
we found a clear dependence on the jet cone size R that can be 
correlated to the suppression of inclusive jets. If part of the en- 
ergy of the medium-induced parton shower is additionally lost 
either through subtraction as soft background or through colli- 
sional processes, significantly larger enhancement in asymmet- 
ric jet production can be obtained for radii of moderate size 
R = 0.4 - 0.5. For small radii, such as R - 0.2, practically all 
of the asymmetry is due to radiative processes. Finally, the full 
calculation has no sensitivity to the choice of the cone radius R. 

The ambiguity of jet/background separation and the diffusion 
of the in-medium parton shower energy through collisional pro- 
cesses 12711 may play an important role in the generation of A j. 
Recovery of the lost energy at large angles outside of the jet 



cone I2O1I would tend to favor the latter scenario. We are cur- 
rently looking for new observables ll40l that have reduced sen- 
sitivity to QGP background fluctuations but retained clear de- 
pendence on the amount of collisional energy loss of the parton 
shower. We finally point out that a statistical peak in the asym- 
metry distribution at finite Ay = 0.4 - 0.5 is likely not related to 
the physics of jet quenching. 

To summarize, significant experimental progress in jet 
physics in heavy ion collisions has been made to date. How- 
ever, taken at face value, the PHENIX and STAR preliminary 
results in Au+Au and Cu+Cu collisions at RHIC Il5lll7ll and 
the ATLAS and CMS jet results in Pb+Pb collisions at the 
LHC ifigl Eoll do not allow, at present, to construct a consistent 
picture of jet modification in heavy ion reactions. Specifically, 
the pronounced dependence of the inclusive jet Raa on the cone 
radius R at RHIC is inconsistent with argued independence of 
the di-jet asymmetry Aj on R. Such independence can easily 
arise from the inherent ambiguity in the separation of the jet 
from the enormous soft background in heavy ion reactions. It 
can also be related to the diffusion of a part of the in-medium 
parton shower energy away from the jet axis due to collisional 
processes. The ATLAS and CMS di-jet measurements have 
played an important role in emphasizing the importance of jet 
studies in heavy ion reactions at the LHC. However, a more 
comprehensive suite of measurements of inclusive jet and di-jet 
modification and further developments in theory are necessary 
to find a consistent picture of jet production and modification 
in dense QCD matter. 
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